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Abstract. Recent A^-body simulations have shown that Einasto radial profiles provide the most 
accurate description of dark matter halos. Predictions based on the traditional NEW functional form 
may fail to describe the structural properties of cosmic objects at the percent level required by preci¬ 
sion cosmology. We computed the systematic errors expected for weak lensing analyses of clusters 
of galaxies if one wrongly models the lens density profile. Even fhough fhe NEW fils of observed 
langenlial shear profiles can be excellenf, viral masses and concenfrafions of very massive halos 
(> 10 ^^Mq//i) can be over- and underesfimafed by ~ 10 per cenf, respectively. Misfiffing effecls 
also sfeepen fhe observed mass-concenfrafion relation, as observed in mulfi-wavelengfh observations 
of galaxy groups and clusters. Based on shear analyses, Einasto and NFW halos can be sef apart 
eifher wifh deep observations of exceptionally massive sfrucfures (> 2 x 10 ^®Mq//i) or by slacking 
fhe shear profiles of Ihousands of group-sized lenses (> 10 ^^Mq//i). 
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1 Introduction 

The hierarchical cold dark matter paradigm with a cosmological constant (ACDM) is highly success¬ 
ful in describing the properties of the Universe and the evolution and formation of structure therein. 
Dark matter halos are traditionally modeled as Navarro-Frenk-White (NFW) density profiles [1, 2]. 
As wifh many ofher ingredienfs of fhe ACDM cosmology, fhe main reasons behind assuming NFW 
halos are fhaf fhey resulf from A'-body simulafions and fhey work relafively well in modeling observ¬ 
able properties of galaxies and clusfers. 

High precision asfrophysics nowadays requires fhaf a level of accuracy of ~ 1-2 per cenf on 
fhe measuremenf of densify profiles of dark maffer halos is achieved. Recenf invesfigafions (see for 
insfance [3-6] and references fherein) have shown fhaf CDM haloes are described besf by fhe Einasfo 
densify profile, which is given by d In p/d In r oc r" wifh a constanf value of a [7] . NFW and Einasfo 
halos differ mainly af very small radial disfances from fhe cenfer of very massive sfrucfures. These 
regions are difficulf fo probe observafionally, and as a consequence if is generally assumed fhaf NFW 
profiles can work as well as Einasto models. 

Problems, however, might arise when highly non-linear processes are involved. Deviations 
of the NFW functional form from A-body results are generally small for halos that are not very 
massive, but can be significant for massive clusters at high redshift [5]. Halos that are well fitted by 
NFW models in 3D may be not NFW-like in projection [6]. This is caused by the halo triaxiality and 
by the effects of substructures and additional matter along the line of sight. Furthermore, observable 
quantities are generally connected to the gravitational potential and its derivates rather than the mass 
distribution. A reasonable characterization of the 3D density profile may fhen fail fo reproduce all 
fhe observable feafures of a halo. A descripfion of fhe halo in terms of the potential could be more 
sensible [8]. 

Well known systematic effects can bias the measurement of mass and concentration in weak 
lensing analyses and the proper comparison to theoretical predictions. The presence of a dominant 
central galaxy has to be accounted for in studies of the cluster core. Baryonic effects can influence 
fhe dark maffer halo profiles [9] . Irregular clusfers may exhibif subsfrucfures and deviafions from a 
simple paramefric profile [10]. Halo concenfrafion can be over-esfimafed in halos elongafed along 
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the line of sight [11, 12]. Cluster member galaxies can be erroneously treated as background source 
galaxies due to wrongly assigned photometric redshifts [13]. Staked analyses can be affected as well 
by scatter in the observable mass proxy and by scatter in concentration at fixed mass. Miscentering 
effect on small scales and the 2-halo term on large scales have to be considered too [14]. Finally, 
selection effects can plague the analysis of not statistical samples of clusters [6, 15]. 

The low level of statistical uncertainties expected in ongoing and upcoming large area surveys 
demands for a proper treatment of all sources of systematic errors. In this paper we quantify the sys¬ 
tematic errors made by modeling the observable weak leasing properties of Einasto halos by means 
of NEW profiles. The NEW models fail mainly in irregular clusters out of equilibrium. Eor these sys¬ 
tems our results will provide a conservative estimate of the minimum error. In our analysis we hence 
assumed that dark matter halos are in fact Einasto-like, and that we have been wrongly modeling 
them with the NEW parameterization. 

Weak leasing analyses provide reliable mass measurements [16-19]. Notwithstanding the long 
history of gravitational leasing, leasing by Einasto haloes is still in its infancy [20, 21, and references 
therein]. The main reason is that the properties of Einasto lenses cannot be expressed in terms of 
simple functions. [22] computed some approximated interpolating functions. [20] and [21] expressed 
the convergence and the shear in terms of Eox H functions. The complementing analysis of the 
leasing properties of the Sersic model can be found in [23]. 

On the observational side, [24] verified that Einasto and NEW modelings gave consistent results 
within errors in the analysis of stacked leasing from galaxies and clusters in the Sloan Digital Sky 
Survey (SDSS). [17] described the stacked tangential shear signal of 20 massive CEASH clusters 
[Cluster Eensing And Supernova survey with Hubble, 25]. They found that the shape parameter, 
mass and concentration of the Einasto model were statistically consistent with the NEW-equivalent 
parameters. The Einasto profile together with other dark matter models was also tested with stacked 
gravitational leasing of high-mass clusters in [26] . 

We adopted as a reference background model the 2013 Planck cosmology [27], with Dm = 

1 — Da = 0.3175, Hubble constant Hq = lOOfi km s“^Mpc“^ with h = 0.671, and us = 0.8344 
for the amplitude of the linear matter power spectrum, which we computed as outlined in [28]. 

2 Halo models 

The Einasto profile [7] provides an excellent fit to dark matter simulated halos, 

/ 2 

PEin = P-2 exp - 

[ a 

where a is the shape parameter, r _2 is the radius where the logarithmic slope of the density profile is 
—2, and p _2 is the density at r_ 2 . The NEW density profile follows the functional form [1], 

PNFW = V / T " T 2 ’ 

(r/rs)(l + r/rs)^ 

where rs(= r_ 2 ) is the scale radius and ps(= 4p_2) is the scale density. 

A halo can be generally described in terms of its mass and concentration. We adopted M 200 , 
the mass inside the sphere of radius r 2 oo> wherein the mean density is 200 times the critical density 
at the halo redshift. The concentration is defined as C200 = F 2 oo/f- 2 - This definition may contrast 
with the common-sense notion that the more concentrated object is the one that has the denser central 
region and the less dense outer halo. In fact, the mass inside the scale radius of the Einasto halo also 
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depends on the shape parameter. An alternative definition is connected to the ratio of the maximum 
circular velocity to the virial velocity, 14iax/f^200 [5, 29]. 

The Einasto and NEW profiles are quanfifafively differenl af small and large radii. On scales 
of inferesf for weak gravifafional lensing (0.1 < r/r 2 oo 5 1-0) the NEW profile can only mimic 
an Einasfo profile wifh a ~ 0.2 [4]. Einasto halos have three parameters, one more than NEW. 
This additional degree of freedom is not the reason behind their better performances [5]. In fact, 
the logarithmic slope is related to the peak height, defined as u{M 2 oo, z) = Sc/cr{M 2 oo, z), and if 
increases wifh mass [30]. In fhe following, we fixed fhe shape paramefer fo u as [5], 


a = 0.115+ 0.0165Z/2. 


(2.3) 


3 Einasto versus NEW lensing 

Eef us now compare fhe lensing properfies of NEW halos fo Einasfo lenses wifh fhe same mass and 
concenfrafion, see Eigure 1 . Concenfrafions were associated fo masses following fhe c — M relafion 
for Einasfo halos derived in [4]. Eor fhe normalization, we considered fhe lens af redshifl Zd = 0.3 
and fhe background sources af Zg = 1.0. 

Differences in convergence, i.e. fhe lensing renormalized projecfed surface densify, are more 
pronounced af large radii whereas differences in fhe reduced fangenfial shear profiles show up af 
small radii. 

Eor M 200 = 5 X 10^^Mq//i, fhe NEW convergence is similar fo fhe frue (Einasfo) convergence 
up fo r < 5 Mpc//i and fhen overesfimafes if. On fhe ofher hand, fhe reduced shear profile is 
underestimated by < 5 per cent between 0.2 < r < 1 Mpc/h. At very small (r ~ 0.05 Mpc//i) or 
very large (r ~ 3 Mpc/h) radii, g^- is overestimated by ~ 10 per cent. 

The lensing properties of NEW and Einasto lenses are very similar to each other in a wide radial 
range for a lens with M 200 = IO^^Mq/ h. The NEW reduced shear reproduces the real profile wifhin 
~ 5 per cenf in fhe range 0.2 < r < 3 Mpc/h. Differences are mosf nofable in very massive halos. 
The reduced shear of lenses wifh M 200 = 2 x 10 ^®Mq//i is overestimated by > 10 per cenf in fhe 
cenfral regions (r ~ 0.3 Mpc/ h) and underesfimafed by < 10 per cenf af r ~ 1 Mpc/ h. On fhe ofher 
hand, fhe convergence is well reproduced up fo r ~ 1 Mpc/h. 

Whereas shear profiles mosfly differ af small radii, fhe convergences are unlike af large radii. 
The combined analysis of galaxy tangenfial disforfion and fhe complemenfary lensing magnificalion 
[17] can be fhen effective in disfinguishing Einasfo from NEW lenses. 

4 Method 

In order fo esfimafe fhe sysfemafic errors experienced by wrongly modeling fhe observed lensing 
properfies of a clusfer, we tiffed fhe properfies of fhe Einasfo halos wifh NEW profiles and compared 
fhe resulfs of fhe regression fo fhe frue parameters. 

We characferized fhe simulafed Einasfo halo of mass M 200 wifh a concenfrafion given by fhe 
median c — M relafion in fhe Planck cosmology [4]. The shape paramefer was relafed fo fhe mass 
wifh Eq. (2.3). As a main fesf we considered lensing observafions of fhe reduced fangenfial shear 
profile p_|_. The function in fhis case is 



(4.1) 
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Figure 1. Profiles of convergence (black), tangential shear (blue), and reduced tangential shear (green). The 
full and dashed lines are for the Einasto and the NEW prohles, respectively. The vertical gray lines bracket 
the radial range between 1 ’ and 15 ’. The vertical red line marks r2oo- The lens is at Zd = 0 . 3 , the sources at 
Zs = 1 . 0 . Top panel: M200 = 5 x IO^'^Mq/H and C200 = 4 . 0 . Middle panel: M200 = 10 x IO^^Mq/Zi and 
C200 = 3 . 7 . Bottom panel: M200 = 20 x IO^'^Mq/Ii and C200 = 3 . 4 . 


where g^- is measured (with statistical uncertainty 5+) in circular annuli at angular position Ex¬ 
pressions for the shear induced by a NEW halo can be found in [31]. Eensing by Einasto deflectors 
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is discussed in [20, 21]. 

We considered observational weak lensing conditions obtainable with deep lensing ground- 
based programs such as CLASH or the Sloan Giant Arcs Survey [SGAS, 32, 33]. We calculated the 
shear in 10 discrete radial bins spanning the range [ 02 D,min = ^ 2 D,max = 15'] with a constant 
logarithmic angular spacing. The redshift of clusters was selected to be Zd = 0.3. The observational 
uncertainty 5+ was computed by assuming a background source population with surface density of 
rig = 20 galaxies per square arcminute at Zg = 1.0 and a dispersion in galaxy intrinsic ellipticities 
of (Te = 0.3. We added the cosmic noise due to uncorrelated large scale structure projected along the 
line of sight following [34]. Non-linear evolution was computed as outlined in [35]. We associated 
errors to each shear measurement but we did not scatter the shear estimates. Thus the expected 
for the true model is zero. 

As a second test, we considered the 3D fit of the density profile. We minimized the function: 

xId = X/ - lnpNFw(G; 4430,200, C3D,20 o)]^ , (4-2) 

i 

where the density p was measured at 50 radial positions from the cluster center spanning the range 
[F 3 D,min = 0.01r2oo, F 3 D,max = f 20 o] with a Constant logarithmic radial spacing. The 3D profile 
is nof an observable like fhe shear profile. Filling procedures like Eq. (4.2) are however roulinely 
employed in sludies of A^-body simulafed samples lo sludy halo concenfralions [4, 6]. 

5 Results 

5.1 Mass and concentration estimates 

The 3D fil consislenlly underestimates fhe masses, see upper panel of Figure 2. The bias is small 
(less lhan 2 per cenf) up fo M200 = 5 x IO^^Mq/ h, bul can be significanl (up lo 10 per cenl) for very 
massive halos (M 200 ~ 2 x 10^^Mq//i). Mass estimates from fhe 3D or fhe 2D fil are very similar 
for low masses (1 < M 200 ^ 3 x 10^^Mq//i). On fhe olher hand, fhe 2D fil overestimates masses 
larger lhan M 200 ~ 5 x 10 ^‘^MQ/h. 

Errors in fhe deferminalion of fhe concenlralion are larger, see bolfom panel of Figure 2. The 
3D fil severely overeslimales C 200 by > 10 per cenf for low mass halos (M 200 rsj IO^^Mq/Zi) and 
significanlly undereslimafes — by nearly 30 per cenl — fhe concenfralions of very massive objecls 
(M 200 ~ 2 X lO^^Mo/Zi). 

Errors expecled in fhe 2D fil of fhe concenlralion are usually smaller bul Ihey can be still sig¬ 
nificant The error is negligible al M 200 10^^ or ~ IO^^Mq/Zi. The maximum overestimate 

(~ 7 per cenl) is for M 200 ~ 3 — 4 x 10^‘^Mq/Zi, whereas fhe concenlralion of clusters wilh 
44200 ^ 5 X 10 ^^Mq/Zi is underestimated. Errors in masses and concenfralions are mosl significanl 
al Ihe high mass fail of Ihe cluster dislribulion (M 200 ^ IO^^Mq/Zi). Very massive clusters are Ihe 
preferential largels of dedicated weak lensing programs. Nearly half Ihe clusters in eilher Ihe CEASH 
or Ihe WlG [Weighing Ihe Gianls, 36] programs have a mass in excess of M 200 rsj IO^^Mq/Zi. 

5.2 Model identification 

Einaslo and NWF models are difficull lo distinguish based on weak lensing properties only. In fad, 
NEW models provide excellenl fils lo Einaslo shear profiles allhough wilh biased masses and con- 
cenlralions. 

We can compute Ihe minimum number of clusters of given mass lhal have lo be slacked lo reach 
positive evidence for model identification, i.e. a difference of 2 for Ihe Bayesian information criterion 
[37], see Figure 3. 
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Figure 2 . Relative systematic errors as a function of mass made when fitting Einasto halos with a NEW 
functional form. The blue and green lines refer to the fits of the reduced shear profile and to the fit of the 3 D 
density profile, respectively. Top panel: Relative errors in mass. Bottom panel: Relative errors in concentration. 


We considered ground- (see Section 4) or space-based surveys. A Euclid-like survey' can reach 
Ue = 0.25 and an effective number density of rig ~ 30 galaxies per square arcminute. For high-z 
lenses at = 1.0, we considered 10 angular bins between 0.5’ and 5’ and a median source redshift 
of 2 ;s ~ 1.5. The inner radius follows the choice made in the analysis of space-based weak lensing 
profiles in [38], whereas the outer radius accounts for the extended coverage obtainable in large 
surveys. 

Even for deep observations, the peculiar lensing signature can be identified in single clusters 
only if the halo is exceptionally massive. Halo identification is difficult at high redshifts too. Com¬ 
pared to the intermediate redshift lenses, the number of high-z lenses to stack is smaller at low masses 
and it is larger for very massive haloes. 


'http://sci.esa.int/euclid/ 
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Figure 3. Number Nc\ of clusters of mass M200 that have to be stacked to distinguish Einasto from NEW shear 
profiles for different surveys and different lens redshifts Zd- 


Our estimate of the number of clusters to be stacked is conservative. The impact of cosmic 
noise, which is a major source of uncertainty at large radii and high redshifts [39], can be reduced by 
modeling the uncorrelated large-scale structure with data readily available in large photometric and 
spectroscopic surveys [40]. Furthermore, stacking increases the signal-to-noise ratio in the very inner 
regions (0 < 1'), which we did not consider in the fitting of low redshift clusters, where differences 
between Einasto and NEW are more prominent. An optimized radial binning more refined in the 
inner regions can also weight more the differences. 

On the other hand, these gains are counterbalanced by off-centering effects and intrinsic scatter 
in the observable proxy of the mass used to bin the lenses, which smooths out the signal [14]. 


5.3 Mass-concentration relation 

Erroneously analyzing Einasto shear profiles with the NEW functional form steepens the observed 
c — M relation, see Eigure 4. The measured relation is overconcentrated at low masses and under¬ 
concentrated at the large mass tail with respect to the input relation. Errors are within the statistical 
uncertainties but the effect is systematic. The bias in the slope estimation is strongly sensitive to the 
inner minimum radius considered in the fitting procedure. We conservatively considered 0min = ■ 

The smaller, 0min> the larger the bias. The steepening is present also for the 3D fit of the density 
profiles. 

The presence of the steepening does not depend on the assumed effective density of background 
sources and the related observational uncertainty on the measured shear. Even in case of infinitely 
accurate and precise measurements, when the estimated concentration (thin blue line in Eigure 4) 
tends to the best fit value (thick blue line), the measured c — M is steeper. 

The effect we measured is larger than the difference between NEW (green-dashed line in Eig¬ 
ure 4) and Einasto (black-full) 3D fits to W-body simulated haloes [4]. This is due to the details of 
the fitting procedure and to the fact that W-body simulated dark matter halos are not exactly Einasto. 
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Figure 4. The mass-concentration relation. The black (dashed green) line is the input (approximated) c-M 
relation derived in [ 4 ] by fitting simulated halos with Einasto (NEW) profiles. The green solid line is the c-M 
relation obtained by fitting the 3 D density profile of the Einasto halo with a NEW profile. The blue lines are 
the c-M relations obtained by fitting the reduced shear profile of the Einasto halo with the NEW model. The 
thick and thin blue lines denote the best fit results (i.e. which approximates the relation estimated in absence 
of observational etTors) and the bi-weight estimator of the marginalized posterior probability distribution. The 
shaded blue area includes the l-tr confidence region. 


6 Conclusions 


Profiles which excellently fit the 3D density distribution of clusters of galaxies can fail to predict the 
observable projected properties of the corresponding halos at the level of 1-2 per cent required by 
precision cosmology. We estimated the systematic errors in the measurement of mass and concen¬ 
tration made by fitting Einasto-like shear profiles with NEW models. Effects can be prominent for 
very massive halos with virial masses in excess of M 200 ~ IO^^Mq/Zt. Misfitting effects can cause 
overestimations of the mass by ~ 10 per cent and underestimation of the concentration by ~ 10 per 
cent. 

Einasto and NEW halos cannot be yet distinguished by present stacked analyses. Accidentally, 
[24] and [17] considered mass scales (M 200 ^ 10^^ and ~ IO^^Mq/Zi, respectively) where the 
difference is minimal. Euture Euclid-like surveys are going to provide the required accuracy. 

The accurate determination of the lens profile can also benefit by combining shear and magni¬ 
fication analyses [17]. Count-depletion measurements using flux-limited samples can constrain the 
convergence profile and help to break degeneracies. 

Misfitting effects can help to solve the problem of observed massive lensing clusters following 
a steep scaling in tension with predictions from the concordance ACDM paradigm [15, 41]. This 
tension can be mostly solved by considering the strong anti-correlation between lensing measured 
mass and concentration [15, 42], the adiabatic contraction of the halos and the presence of a dominant 
brightest cluster galaxy [9, 15, 43, 44], and, mostly, selection effects [6, 15, 45]. The effect we 
considered goes in the same direction. 
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